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Abstract 

In this work we analyze the stochastic dynamics of the Kauffman 
model evolving under the influence of noise. By considering the av- 
erage crossing time between two distinct trajectories, we show that 
different Kauffman models exhibit a similar kind of behavior, even 
when the structure of their basins of attraction is quite different. This 
can be considered as a robust property of these models. We present 
numerical results for the full range of noise level and obtain approxi- 
mate analytic expressions for the above crossing time as a function of 
the noise in the limit cases of small and large noise levels. 
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1 Introduction 



The Kauffman model (or N-K model) describes the dynamics of a network of 
N Boolean spins, each controlled by K other spins through a binary function. 
It was first proposed by S. A. Kauffman in 1969 |lTJ as a model for cell 
differentiation and genetic networks. Since then, its application has been 
extended to many other fields in physics, biology, computational and social 
sciences. During the past few decades, most of the work done on Kauffman 
models has been dedicated to the study of the configuration space structure, 
the length and number of cycles, the size of basins of attraction, and the phase 
transition between ordered and disordered phases (for references see ]![). 
These properties are obtained by considering the deterministic dynamics of 
the system, which is well known by now. However, those studies have shown 
that some of the generic properties of the Kauffman model are far from being 
robust. The non- robustness of the deterministic dynamics is reflected, for 
example, in the fact that by slightly changing a given initial configuration of 
spins, the system may "jump" from one basin of attraction to a very different 
one. On the other hand, due to the exponential growth of the state space 
with N, it is often necessary to thoroughly probe the state space in order 
to determine a generic property of the system, such as the mean number 
of different basins of attraction or the mean cycle length. Actually, it has 
recently been shown that a systematic bias due to an under-sampling of the 
state space can be present in some of the results reported in the literature 
during the last 30 years || . Therefore it is valuable to find a method which 
reveals the robust properties of the Kauffman model. 

Real networks are always subjected to external fluctuations. Conse- 
quently, the relevant properties characterizing the network should exhibit 
a certain degree of robustness to external perturbations. In 1989, both Mi- 
randa et al. EI and Golinelli et al. H analyzed the stochastic dynamics of 



the Kauffman model in the case in which an external noisy signal is present. 
In this work we extend the study of the stochastic dynamics of the Kauffman 
network with noise, by performing more accurate numerical simulations as 
well as analytic calculations. We focus our attention on the time it takes for 
two trajectories, starting out from different initial conditions, to cross. We 
consider two cases. First the situation in which each one of the N spins is 
determined by K other spins chosen randomly from everywhere in the sys- 
tem (the Kauffman net). The second case is a d-dimensional lattice in which 
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each spin is preferentially coupled to its immediate neighbors. As we will 
see, both models exhibit qualitatively similar behavior. 

In section 0, we introduce the Kauffman model with deterministic and 
stochastic dynamics. In section |3] we describe our numerical results for both 
Kauffman nets and lattices. In section f|, we study closely the behavior of 
these models in the limits of small and large noise. We summarize this work 
in section |5] with a brief discussion of the results. 

2 The Kauffman model 

2.1 Deterministic dynamics 

A Kauffman model consists of iV Boolean spins {Si, S 2 , ■ ■ ■ , S^} with Si being 
either zero or one. The value of each spin Si at time t + 1 is determined by 
the values of K other spins S^, Si 2 , . . . , Si K , which are called the controlling 
elements for spin Si. (The number K is called the connectivity of the system.) 
Once the connections in the system are established, each spin Si is assigned 
with a Boolean function fa of its K controlling elements. A realization of 
the Kauffman model consists of the set of connections and Boolean functions 
assigned to every spin. The dynamics of the network is then given by 

Si{t + 1) = fi(S h (t), S l2 (t), S iK (t)) for i = 1, . . . , N. (1) 

For convenience, we will denote by Et the state of the system at time t: 

E t = {S 1 (t),S 2 (t),...,S N (t)}. 

In different Kauffman models, the assignments of the K controlling spins 
S^, Si 2 , . . . , Si K of each spin Si and the dynamic rules are different. In 
Kauffman nets, the controlling elements of Si are assigned randomly, whereas 
in a Kauffman lattice they are chosen only among its nearest neighbors. The 
dynamic rules /j are chosen randomly in such a way that its two possible 
outcomes, and 1, occur with probability p and 1 — p respectively. If the re- 
alization of the network is time-independent, the network is called quenched, 
while if either the set of connections or the set of Boolean functions fi are 
re-assigned at every time step, the network is termed annealed. 

Annealed models are more convenient for theoretical studies than quenched 
models. For example, by using the annealed approximation it has been shown 
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analytically that Kauffman nets exhibit three different phases: frozen, criti- 
cal and chaotic, depending upon the values of the parameters K and p [§]. 
The critical value of the connectivity is given by K c = [2p(l — For 
K < K c the system is in the frozen phase, whereas if K > K c it is in the 
chaotic phase. Throughout this work we will use p = 1/2, for which K c = 2. 

But for most real cases (neural networks, genetic networks, etc.), quenched 
models will be more appropriate since in real networks neither the connec- 
tions nor the interactions between the elements change randomly at every 
moment. However, it has been shown that in the limit N — > oo, both the 
quenched and the annealed Kauffman nets are exactly equivalent with respect 
to the evolution of the overlap between different configurations, although not 
with respect to the configurations themselves [TO], M. In this paper, our 



main focus is on quenched Kauffman nets and lattices. 

Due to the finite size of the system, there are a finite number of possible 
configurations, to wit Q = 2 N . Therefore, starting out with an initial config- 
uration, the system will eventually fall into a previously visited state, after 
which the same sequence of states repeatedly occurs again. The state space 
breaks up into a multitude of cycles (or attractors). The totality of points 
which end up in the same attractor represents its basin of attraction. 

2.2 Stochastic dynamics 

There are different ways of introducing noise into Kauffman models and they 
reveal different features of the configuration space of the model. Following 



Miranda and Parga [12], we introduced noise in the following way: 



s ./j +1 \ = ( /"(^iW'^Wv.^W) with probability 1 - r, 
\ 1 - fiiSi^t)^^), . . . ,S iK (t)) with probability r. 

In this way, every spin S{ has a probability r of violating the deterministic 
rule (p. We will say that an n-spin flip event has occurred at a particular 
time step, if n spins violated the deterministic rule in this time step. Notice 
that this stochastic dynamic rule has a symmetry about the point r = 0.5. 
For r > 0.5, if we make a substitution fi — ► 1 — fi, the rule becomes identical 
with the case of 1 — r. Since the Boolean functions fi are assigned randomly, 
fi and 1 — fi are equally likely to appear in a particular realization of the 
model. After averaging over different realizations, the cases with probabilities 
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r and 1 — r are indeed identical. Due to this symmetry in the stochastic 
dynamical rule (H), we only need to consider the case r e [0, 0.5]. Note that 
for the particular value r = 0.5, Si(t + 1) is equally likely to be zero or one 
independently of the value of /j. 

In the presence of noise, the concept of "attractor" does not hold any 
more, since as the system evolves, there is a non-zero probability of "jumping" 
to a different attractor, and consequently every point in the state space can be 
reached from any initial condition. In this sense, the "boundaries" between 
different attractors become more diffuse as the level of noise increases 
However, we shall argue that the system has a sort of effective attractor even 
in the presence of noise, specifically when r is large. 

One of the interesting things to study is the time it takes for two trajecto- 
ries to cross. Suppose that we start with two different initial configurations, 
E and E , and let them evolve according to (0), noting all the configurations 
produced: 




The crossing time r is then defined as the time for which either one of the 
trajectories coincides for the first time with a configuration previously vis- 
ited by the other trajectory. For example, when E r is equal to any of the 
configurations E , E l5 E 2 , . . . , E T . Two important cases have to be distin- 
guished, when E and E belong to the same basin of attraction, and when 
they belong to different basins of attraction. For those cases we will denote 
the crossing time by r s and r^, respectively. 

Miranda and Parga examined the behavior of the system by considering 
only the attractors with largest and next-largest basins of attraction. They 
then showed that for small values of r, the behavior of r s and are very 
different. At r = 0, two trajectories from the same basin will cross in a time 
comparable with the sum of two times: first, the transient time required to 
enter the attractor, and second, the length of the attractor itself. Conversely, 
two trajectories starting out from different basins will never cross. On the 
other hand, they found that for sufficiently large values of r, the crossing 
time became independent of the starting point. It did not matter where the 
two trajectories start, the two basins merge into a sort of effective attractor 
and the trajectory bounces around within that subset of the system-states. 
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For these larger values of r, the observed effective basin size increased with 
r. From this, they drew the conclusion that the disappearance of basins of 
attraction with the increase of r is a sort of hierarchical process in the sense 
that in a finite period of time, the portion of the whole state space explored 
by a trajectory starting out from a given basin of attraction increases with r. 
Complete randomness is achieved at r = 0.5, where the trajectory explores 
the entire state space. 

As we will see, our simulation will show the same general behavior as 
described by Miranda and Parga. But we shall explore the behavior in more 
detail, showing the crossing time for the whole range of values of r and K. 

3 Numerical results 

Kauffman nets and Kauffman lattices differ in the structure of their basins of 
attraction. One would expect this difference to be reflected in the response 
of these models to the influence of noise. For random realizations of the 
coupling functions /«, what determines the basin structure is the connectivity 
K. Therefore, we will first analyze separately the cases with large K (chaotic 
phase) and small K (ordered phase). 

We will partially follow Miranda and Parga's approach in that we compute 
the average crossing time t s by using two initial configurations, S and S , in 
the largest basin of attraction. For the average crossing time r^, we pick one 
starting configuration in the largest basin and the other in the next largest 
one. The reasons to choose only the two largest basins of attraction will be 
clear in what follows. 

3.1 Kauffman models with large K 

We want first to characterize the structure of the basins of attraction. One 
way of doing it is by computing the distribution of basin-sizes W(n), which 
is the fraction of the state space Q occupied by the n-th largest basin of 
attraction. In Fig. p] we show W(n) for a Kauffman net and a 1-dimensional 
Kauffman lattice, both with N = 20 and K — 5. It can be seen from this 
figure that both models exhibit a very similar structure in their basins of 
attraction in the sense that the basin sizes are similar. It is worth mention- 
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ing that in other aspects, like orbit length or transient timeQ, the basins of 
attraction can still be very different in both models. 

From Fig. |I| it also can be seen that the largest and next largest basins 
occupy more that 90% of the whole state space. Therefore, to a good ap- 
proximation it can be assumed that the dynamics takes place mainly in these 
two largest basins. 
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Figure 1: Distribution W(n) of basin-sizes for: (a) the one-dimensional Kauffman 
lattice and (b) the Kauffman net, both with N = 20 and K = 5. For the Kauffman 
net, the connections between spins are chosen randomly, whereas for the lattice 
every spin is connected to itself and to its 4 nearest neighbors (periodic boundary 
conditions were used). The number n in the horizontal axis corresponds to the 
n-th largest basin in the model. 

Fig. shows the average crossing times t s and as functions of r for the 
Kauffman net and the 1-dimensional lattice both with N = 20 and K = 5 
(chaotic phase). Notice that these two kinds of N-K models exhibit very 
similar behavior under the influence of noise. 



The transient time is the time it takes before a trajectory enters the stable cycle. 
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Figure 2: Average crossing time t for different Kauffman models with N = 20 and 
connectivity K = 5 (chaotic phase) as a function of (a) the noise intensity r and 
(b) the inverse of the noise intensity. The symbols are as follows. 1-dimensional 
Kauffman lattice: (□) and (o) t s . Kauffman net: (+) and (*) T s . Each point 
is the average over 4000 realizations of the model. 



3.2 Kauffman models with small K 

The ordered phase is characterized by K = 1 and K = 2. In this section 
we will present the results for the minimum value of K, namely, K — 1. 
In Fig. m we show W(n) for a Kauffman net and a 1-dimensional Kauffman 
lattice, both with N = 20 and K = 1. The connections in the Kauffman 
net were, as usual, chosen randomly, whereas in the 1-dimensional lattice the 
node Si was connected either to SV-1 or to Si + i with equal probability (we 
use periodic boundary conditions). 

From Fig. |3] it is apparent that in this case, the basin structures of the 
Kauffman net and the Kauffman lattice are less similar than in the chaotic 
phase. For the lattice, the two largest basins of attraction no longer occupy 
more than 90% of the whole state space, whereas in the Kauffman net they 



8 



still do. However, the response to the influence of noise is mostly the same 
in both models, as can be seen from Fig. f| where the crossing times t s and 
Td are plotted as functions of the noise intensity r. 
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Figure 3: Distribution W(n) of basin-sizes for: (a) the one-dimensional Kauffman 
lattice and (b) the Kauffman net, both with N = 20 and K = l(frozen phase). 
The number n in the horizontal axis has the same meaning as in Fig. |l|. 



3.3 Robust behavior of the crossing time 

Figures ^| and [| show that the behavior of the different Kauffman models 
under the influence of noise have the following general characteristics, both 
in the frozen and in the chaotic phases: 

• For small r, decreases as 1/r while t s is nearly constant. 

• For large r, both Td and t s increase with r and become equal at r = 0.5. 

• For intermediate values of r, has a minimum when t s w Td- 

Let us analyze separately each one of the above characteristics. 
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Figure 4: Average crossing time r for different Kauffman models with N = 20 and 
connectivity K = 1 (frozen phase) as a function of (a) the noise intensity r and 
(b) the inverse of the noise intensity. The symbols are as follows. 1-dimensional 
Kauffman lattice: (□) and (o) t s . Kauffman net: (+) and (*) t s . Each point 
is the average over 4000 realizations of the model. 



In the limit r — > 0, r s approaches a finite value t s (0), the mean crossing 
time for two trajectories in the largest basin of attraction in the absence 
of noise. This crossing time is roughly one half the average cycle length, 
plus one half the average transient time. The transient time and the cycle 
length are of the same order of magnitude, therefore r s (0) is expected to 
be approximately equal to the mean cycle length (L) of the largest basin of 
attraction (see Fig. ^a). 

On the other hand, diverges as r — > 0. The numerical data (see fig- 
ures 0b and £§b) suggest that in this limit, has the form 

Td ^a{K 1 N)_ + b ^ N ^ r _ >()> (3) 
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The above divergence is due to the fact that, in the absence of noise, there 
is a zero probability for a trajectory to jump between different attractors. 
Under the deterministic dynamics, every trajectory will remain within its 
own basin of attraction for ever. In the next section we will see that the 1/r 
behavior of is a consequence of the fact that the dynamics is governed by 
one-spin flip events when r is small. For large values of K, the largest basin 
occupies almost the whole state space. Under these circumstances, every 
time a one-spin flip occurs, the trajectory in the next largest basin will have 
a finite probability of diverging very substantially from the path it would 
have followed in the absence of noise. That divergence will usually force the 
trajectory into the largest basin. In fact, for a fraction of order one of the 
noise events in the smaller basin, the noise will flip the trajectory into the 
largest one. Once the two trajectories are in the largest basin, they have a 
lifetime b(K, N) before they cross. This lifetime is expected to be of order 
one of t<,(0), the typical length of the largest basin of attraction. The above 
can actually be seen in Fig. || from which it is apparent that for large K, 
b(K, N) « r s (0) « (L). 

In the opposite limit r — > 0.5, the crossing times r s and Td become equal, 
which means that for high levels of noise, the barriers between different at- 
tractors become small. When r reaches its maximum value 0.5, all barriers 
vanish. In this case, both trajectories randomly jump from one state to an- 
other throughout the state space, and both t s and become equal to the 
time it takes for two random walks to cross. As derived in section |4.2| , this 
crossing time is the solution to the "birthday problem" , i.e. 

r s = r d oc 2 N ' 2 (4) 

In this way we have obtained a qualitative description of the limiting 
cases of figures |2| and £|. The one qualitative feature left to describe is the 
crossover from the small r to the large r behavior. As one can see from these 
figures, the crossover occurs when the two times r s and r d become roughly 
equal. This in turn happens when 

r ~ a(K,N)/r s (0) ~ 1/JV (5) 

Thus, the minimum in occurring between the two previous limit values 
of r, can be interpreted as the result of a sort of "competition" between the 
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Figure 5: (a) Plot of the mean cycle length (L) as a function of r s (0) for a 
Kauffman net with N = 12. Each point corresponds to a different value of K, 
starting with K = 1 for the first point in the lower left corner of the graph and 
ending with K = 12 for the last point in the upper right corner. The dashed line is 
the best linear fit to the numerical data (circles), (b) Same type of graph as before 
but now showing the dependence of b(K,N) on r s (0). The slopes of the dashed 
lines are (a) ~ 0.96 and (b) ~ 0.90, which shows that b(K, N) w r s (0) w (L). 



randomness in the system (coming from the presence of noise), and the bar- 
riers separating the attractors (which come from the deterministic dynamics 
of the system). 

3.4 A Kauffman model with equal basin-sizes 

Finally, we would like to mention that the above results are also true for 
Kauffman models in which all the basins of attraction have the same weight. 
As an example, consider a Kauffman lattice with N = 20 and K — 1 in which 
every spin is connected to itself. For K = 1 there are only four Boolean func- 
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tions fi(S): tautology fi(S) = 1, contradiction fi(S) = 0, identity fi(S) = S 
and negation fi(S) = 1 — S. Imagine then the very specific realization in 
which two of the coupling functions are identity two are negation and all 
the others are either tautology or contradiction. By simple analysis, we 
know that for this specific model, the whole state space is composed of eight 
basins of attraction with equal size. Fig. ^ shows the crossing time r as a 
function of r for this particular model. Since in this case all the basins of 
attraction have the same size, the two initial conditions needed to compute 
r were chosen randomly among the whole state space. Again, the r ~ 1 jr 
behavior for small r and the r oc 2 N / 2 behavior for r — > 0.5 are obtained. Of 
course we can construct many other models by choosing different coupling 
functions fi for this K = 1 self-correlated case. All the Kauffman models we 
have explored have shown this kind of behavior under noise. 

4 Theoretical analysis 
4.1 Small r limit 

In this region, the main characteristic of average crossing time is that ~ 
1/r. The reason for this dependence is that the stochastic dynamics is dom- 
inated by one-spin flip events, in the following sense. The probability of a 
one-spin flip event (~ r) is much larger than the probability of a two-spin 
flip event (~ r 2 ). If the probability to jump to a different basin of attraction 
in a one-spin flip event is significantly different from zero, then the dynamics 
will be dominated only by this kind of events. Even if it was necessary to flip 
two spins to jump from one basin to another, this process can be decomposed 
into two one-spin flip events occurring sequentially, instead of being carried 
out at once in one two-spin flip event. 

As we have shown in section || the 1/r behavior is present in a wide 
variety of Kauffman models. This in turn, implies that the one-spin flip 
events dominate the dynamics for small values of the noise. In this subsection, 
our goal is to derive an expression for the coefficient a(K, N) for Kauffman 
nets. To do so, we will make the assumption that the dynamics takes place 
only in the two largest basins of attraction. Although the ~ 1/r behavior 
is generally true, as we have found, the preceding assumption is not true for 
all Kauffman models, especially for those with small values of K, but as we 
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Figure 6: Average crossing time t(O) for a Kauffman lattice with TV = 20 and 
connectivity K = 1 as a function of (a) the noise intensity r and (b) the inverse 
of the noise intensity. For this model, every spin is correlated to itself and we 
choose the Boolean functions in such a way that the whole state space is composed 
of 8 basins of attraction with equal size. Each point is the average over 10000 
realizations of the model. 

show below, it becomes more valid as K increases. 

To start the calculation of a(K, N), let P\ 2 be the probability for jumping 
from the largest basin to the next largest basin with a one-spin flip event 
and P2,i be a similar probability but jumping in the opposite direction^. Let 
us also define Q\ and Q2 as the probabilities of remaining in the largest 
basin and in the next largest one, respectively, after one-spin flip event. 
Simulations show that for Kauffman nets these probabilities have a slight 
dependence on N but a very strong dependence on K. The result is that 
P\,2 + Q\ ~ -^2,1 + Q2 ~ M, where M is approximately constant for all values 

2 By definition, the ratio of P\ t 2 and P?^ is strictly the inverse of the ratio of the size 
of the largest basin to the next largest one in one realization. 
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of K and M > 0.92 (see Fig. 0). 

However, the dynamics depends not only on the total sum M but also 
on the particular values of Pip, Qi, Pi,\ an d Qi- For small values of K, 
both P 12 and P 2j i are very small and comparable with (1 — M). So, even 
though M is large, the interaction between the two largest basins is weak in 
the sense that the probability of jumping into smaller basins is of the same 
order as Pip and P%i- Therefore, for small values of K the smaller basins 
play a significant role in the dynamics of the system. The above can be 
seen in Fig. [7|, in which the probabilities Pip, P 2) i, Qi, Q2, and the sum M 
are plotted as functions of K. From this figure it is apparent that the two 
largest basins of attraction are the dominant ones for large values of K (say 
K > 5). The closeness of M to 1 means that a trajectory will seldom jump 
into a basin other than the two largest ones with only one-spin flip event. In 
view of this result, in some of the arguments below we will assume that K 
is sufficiently large so that the dynamics takes place only in the two largest 
basins of attraction. 

With the information about these probabilities, we can give an approx- 
imate calculation of t& as a function of r for the case in which K is large. 
We know that every spin violates the deterministic rule (|l|) with probability 
r. Therefore, the probability of a one-spin flip event is Nr and consequently 
the expected time for this event to occur is T = l/(JVr). This is true for 
all configurations. For sufficiently small values of r, this expectation time 
is much longer than the average crossing time for two configurations in the 
same basin. Hence, once two configurations jump into the same basin, their 
trajectories meet before the next spin-flip event becomes possible. 

There are two cases in which the two trajectories meet after the occur- 
rence of a one-spin flip event at time T: the configuration in the largest basin 
remains in it while the configuration in the next largest basin jumps into the 
largest one, or vice versa. The above occurs with probability {QiP 2 ,i + 
QzPip)- Similarly we get the probability for the crossing of the two trajec- 
tories after one-spin flip events at 2T, 3T, etc. The average value of r d is 
then: 

r d « T • (QxPa,! + Q 2 Pip) + 2T ■ {QiQ 2 + PipP 2 ,i) • (QiP 2 ,i + ^1,2) + • • • 

00 

= ]T m ■ T ■ (Q!P 2il + Q 2 Pip) ■ (QiQ 2 + PipP^T" 1 

m=l 
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Figure 7: Plot of the probabilities P 1)2 (O), P 2 ,i (A), Q 1 (□) and Q 2 (x) as 
functions of K, for a Kauffman net with A?" = 12. Also shown is M as a function 
of K, obtained as M = Pi 2 + Qi (*), and as M = P 2 1 + Q 2 (°). Note that even 
though M is a constant for all values of K, the probabilities Pi 2 and P2,i are 
rather small for small values of K. 



Nr ' (1 - Q1Q2 - Pi, 2J P2,i) 2 
q(JC, AO 
r 

where a (if, AT) is explicitly given by 

1 Q\Pz,X + Q 2 Pl,2 



(6) 



a(if, N) 



N(l- Q!Q 2 - Pi, 2 P 2 , 



(7) 



In the above derivation, which is true for the case in which K is large, we 
have ignored the time b(K, N) it takes for two trajectories to cross after 
they have jumped into the same basin. Fig. |8] shows the coefficient a(K, N) 
obtained from equation (0) and from simulation, for different values of K in 
a Kauffman net with A^ = 12. It can be seen that the simulation and the 



16 



theoretical result agree very well for K > 5. It is worth emphasizing that for 
small values of K and other Kauffman models where the two largest basins 
don't have such dominance, the effect of other basins besides the largest and 
next largest ones has to be considered to perform an accurate derivation for 
the coefficient a(K,N). 

0.5 i . . . . . . . . . . . . 1 
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Figure 8: The coefficient a(K, N) in Eqn. (0) for different values of K for 
a Kauffman net with N — 12. (x) is the result obtained from simulations 
by sampling the whole state space. (□) is the result obtained from the 
theoretical result Eqn. (|7]). 



4.2 Large r limit 

When r acquires its maximum value 0.5, the barriers between different at- 
tractors vanish and the dynamics transforms into a random mapping of the 
state space into itself || f|. In this limit, the crossing times r s and r d be- 
come indistinguishable and we will refer to both of them simply as r. To 
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understand the limit r — > 0.5 we will first give a simple "birthday-problem" 
argument to obtain the order of magnitude of r. Then we will proceed to 
a more elaborate analysis to obtain an approximate expression for r as a 
function of r, valid for r close to 0.5. 

Imagine two walkers moving at random through a space of dimension 
Q = 2 N . They go through a number of steps n — 1, 2, 3 . . . As each step is 
completed, they have a total number of chances C n on landing on some place 
previously covered by the other walker. At the first step we could have the 
two walkers at identical positions. Thus, C\ = 1. At the second step, C2 = 
3 + C\ since each can land on the original position of the other or they can 
both land at precisely the same place. After n steps C n = + 2n + 1 = n 2 . 
For large fl and small n, the probability of having collided with the path is 
then of order p n « C n /Q. The average time for crossing is then roughly given 
by the n-value for which p n becomes of order unity, so that n 2 is of order Q 
or n = 0(tt 1 ' 2 ). Therefore, the crossing; time r satisfies r oc Vt ' . 

To derive a more precise functional relation between r and r we have 
to compute the probability p c (t) for two trajectories {S , Si, . . . , S T } and 
{So, Si, . . . , S T } to cross at time r. In order to do that, we have to have 
first the probability p for two configurations S 4 and Sj to be the same (note 
that t and t might be different). Let Si(t — 1) be in St_i and Si(t — 1) be 
the corresponding spin in Since these spins are in the same position 

(each in its respective configuration), the deterministic rule fi they obey is 
the same. Notice that the dynamical equation (§) can be written as 



From this expression it follows that the probability for Si(t — 1) and Si(t — 1) 
to evolve according to the deterministic rule fi is (1 — 2r) 2 . Let us denote 
by pi the probability that Si(t) = Si(t) when both spins are updated ac- 
cording to the deterministic rule f\. To calculate p\ we follow the annealed 
approximation introduced by Derrida and Pomeau , which leads us to the 
following two possibilities: 

1. The K inputs of Si(t — 1) and Si(t — 1) are the same, which occurs with 
probability 1/2 K . When this happens, Si(t) = Si(t) with probability 1. 




fi with probability 1 — 2r, 

evolve randomly with probability 2r. 



(8) 
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2. At least one of the inputs is different, which occurs with probability 
(1 — 1/2 K ). In this case, if the evolution rules fi are assigned in a 
sufficiently random way, there is a probability of 1/2 that = Si(t). 



From the above it follows that 

Pi = (1 - 2r) 2 



2 K 



+ 1 



1 

2~k 



On the other hand, the probability p 2 that Si(t) = Si(t) when the evolution 
rule fi is violated in one or both of the configurations, is simply given by 

p 2 = [1 - (1 - 2r) 2 ] ■ ~ 

Combining the values of pi and p 2 given above, the probability p for both 
configurations E t and E^ to be equal is 



P 



(1 - 2r) 2 



1 1 1 

2^ + ^ ~ 2^ ' 2. 

N 



+ [l"(l-2r) 2 ]-i} 



^•<! 1 + 7TF-( 1 -2r) 



2^ 



(9) 



If q(t) is the probability that the two trajectories have not yet crossed at 
time t, then p c {t), the probability for the two trajectories to cross at time t, 
is given by 

p c (t) = q(t - 1) - q(t) (10) 

The two trajectories are still separated at time t if none of the configurations 
{S , Ei, . . . , E t } is equal to any of the configurations {E , Ei, . . . , E 4 }. The 
probability for this to happen is 

q(t) = (l-pf + V 2 

Substituting this value of q(t) into equation (ffDj) we get 

p c (t) = (l-pf -(l-p)(^) 2 

Therefore, the average crossing time r = Y^Li t • Pc(t) is given by 

Et-[(l-pf-(l-p) (t+1)2 ] 



T 



(11) 



(12) 



t=i 
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d(l-p) 



dt 



dt 



71 



In (1 — p) 



Expanding the logarithm in the above equation around p — 0, and retaining 
only the terms up to the the first order, we finally get 

nz r o i N / 2 



argument for the case r = 0.5. When r is not exactly 0.5, we do not have 
a simple birthday problem because the coupling between different elements, 
and consequently the functions fi, still play a role in the dynamics. However, 
it is clear from the above equation that the problem can be viewed as a 

r „ iN 

birthday problem with an effective state space Q e ff = 1+ ( 1 _ 2r )2/ 2 jf 

Fig. § compares the theoretical result fll3|) with the numerical simulation 
for K = 1 and K = 5. As can be seen, the analytic result approximates very 
well the numerical data in the region of r close to 0.5. In this region the 
annealed approximation holds because the noise breaks the correlations be- 
tween the spins. However, for small values of the noise those correlations are 
important and cannot be neglected. We therefore do not expect agreement 
between the numerical and theoretical results for small values of r since in 
this region the annealed appriximation is not longer valid. 

5 Conclusions 

We have considered the effect of external perturbations (noise) in the dy- 
namics of the Kauffman model. The behavior of both, the Kauffman net 
and the Kauffman lattice under the influence of noise is very similar, even 
though these models might have a quite different structure in their basins of 
attraction. In this sense, the response of the Kauffman models to the effect 
of noise can be considered as a very robust property. 

In the limit r — > 0, the most important property is the 1/r behavior of 
the crossing time r^, which has been always present in the Kauffman models 
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Figure 9: Average crossing time r as a function of r for two Kauffman nets 
with N = 20 and connectivities (a) K = 1 and (b) K = 5. In both graphs 
the dashed curve is the result of the numerical simulation for t^, whereas the 
solid line is the plot of the analytic expression fllBp. 



we have studied so far. This 1/r behavior is a consequence of the fact that, 
for small values of the noise, the dynamics is dominated by one-spin flip 
events. An approximate equation relating and r was obtained by taking 
into account the fact that, for large values of the connectivity, most of the 
dynamics takes place in the largest and next largest basins of attraction. 

In the second limit r — > 0.5, the barriers between different attractors 
disappear and the dynamics transforms into a random mapping of the state 
space into itself. As a consequence, r s and become equal. In the case in 
which r = 0.5, the crossing time between the two trajectories can be seen as 
the solution of a "birthday problem" in a space of size Q = 2 N . For other 
values of r, but still close to 0.5, the correlations between spins have to be 
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taken into account, which have the effect of reducing the size of the region 
of the state space explored by the dynamics. 

Between these two limit cases for the noise, there is a minimum in the 
value of r d as a function of r. In a loose sense, this minimum could be inter- 
preted as the result of a "competition" between the randomness generated by 
the noise, which tends to homogenize the state space by diminishing the bar- 
riers across different attractors, and the deterministic dynamics, which tends 
to confine two trajectories within the same basin. To analyze this region it 
would be necessary to consider multiple-spin flip events as well as long-time 
step correlations. 

The results and techniques presented in this work could be extended to 
other systems acting under the influence of noise in order to provide them 
with a robust characterization. 
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